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Constraints on dark energy models from radial baryon acoustic 

scale measurements 

Lado Samushia 1 ' 2 and Bharat Ratra 1 
ABSTRACT 



We use the radial baryon acoustic oscillation (BAO) measurements of 
Gaztanaga et al.l (j2008al ) to constrain parameters of dark energy models. These 
constraints are comparable with constraints from other "non-radial" BAO data. 
The radial BAO data are consistent with the time-independent cosmological con- 
stant model but do not rule out time- var ying dark energy. When we combine 
radial BAO and the iKowalski et al.l (120081 ) Union type la supernova data we get 
very tight constraints on dark energy. 



Subject headings: cosmological parameters 
ture of universe — supernovae 



distance scale — large-scale struc- 



1. Introduction 



Recent type la supernova data (IKowalski et al.l 120081 ; iRubin et al.l 120081 ; ISahni et al. 



20081 ) confirm, at high significance, that the cosmological expansion is currently accelerating. 
If we assume that general relativity is valid on cosmological length scales, this cosmological 
acceleration requires that the Universe's current energy density budget is dominated by far 
by an approximately spatially uniform component - dark energy - with negative pressure 
p < — p/3, where p is the dark energy density. 

The most econom ic and the oldest form of dark energy is Einstein's cosmological con- 
stant A (IPeeblea Il984j ) which is time independent and has an equation of state p^ = — Pa- 
The time-independent cosmological constant model (ACDM) provides a reasonably good fit 
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to most current cosmological data (see, e.g. jRatra &: Vogeleyll2008l ; iFrieman et al.ll2008l ). but 
despite this success a lot of alternative models of dark energy have been proposed over the 
years. One reason for this is that cosmological data can not yet tightly constrain the various 
options currently under debate, although it is thought that in next decade a large amount of 
more precise new data will tightly constrain or measure d e partures from the now standard 



ACD M model (see, e.g., IPodariu et al.ll2001al : IWand 120071 : Barnard et al.ll2008l : iTang et al. 



2008hfl The second reason is that a cosmological constant is not straightforward to under- 
stand on a more fundamental level; in particular, it is difficult to accept as fundamental the 
needed new energy scale of a few meV. 

In some scalar field dark energy models a nonlinear attractor solution ensures that this 
scale of a few meV is not fundamental, rather it follows from a much high er energy scale as 



the s c alar field energy densit y decreases during the cosmological expansion (IPeebles fc Ratra 



19881 ; iRatra fc Peebles! 1 19881 ). In this paper we also consider this 0CDM model in which the 
dark energy is modeled as a slowly rolling scalar field <fi with self-interaction potential energy 
density V(4>) oc (f)~ a , where a > 0. A number of other models with time-varying dark energy 
have been proposed^ We also consider the XCDM parametrization which is often used to 
describe time- varying dark energy models^] In XCDM dark energy is assumed to be a perfect 
fluid with effective equation of state p x = u> x p x , where u x is a number less than —1/3. 

Given a model it is possible to compute quantities such as the Hubble parameter H(z) 
and the angular diameter distance cIa(z) as a function of redshift z. Since these quantities 
are observable (or observational results depend on them), given perfect data it is possible in 
principle to compare observational results with theoretical predictions and determine which 
dark energy model provides a better fit to the data. There are a few difficulties however. 
Theoretical predictions not only depend on dark energy model parameters, they also can 
depend on a number of other cosmological parameters, such as the energy density of non- 



1 The ACDM model is based on the cold dark matter (CDM) structure formation m odel which might 
have trouble accounting for some observations (for a summary of some of the problems see IPeebles fe Ratra 
20031) . 



Modifications of general relat ivity on cosmological length scales that might do away with dark energ 



have also been discussed (see, e.g.JWei fc Zhang 



2008 



Tsuiikawa fc Tatekawa 



2008 



Capozziello et al 



200 



Sotiriou fc Faraonl 2008 : Bamba et al. 20081 ). In the models we consider in this pa per, dark energy affect s 



other fields o n ly thr o ugh gravity. For d ar k energy mode l s with other coup li ngs s ee lLa Vacca et al.l ([20091 ) . 



Wu fc Zhan3 (|2008l ). Ijesus et al.l |2008h. Untusch et all ( 20081). IWu et al.l (l2008h. and references t herein . 
For other dark energy models see Izhang fc Chenl (|2008h . TAlcaniz et al.l |2009h . ISetare fc Saridaki] |2009l ). 
Basilakos fc Plionisl ( 2008 ). Shapiro fc Solal ( 2008 ). and references therein. 



3 The XC DM parametrization cannot describe time- varying </>CDM dark energy at late times (see, e.g., 
Ratralll99ll) . 
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relativistic matter or baryonic matter or radiation, the energy density perturbation spectral 
index, the total mass of neutrinos, the Hubble constant, etc. While some of these parameters 
are directly measurable, some of them have to be simultaneously determined from the data. 
This leads to degeneracies in cosmological parameter space. 

The other issue is possible systematic errors in cosmological observations, which are 
very difficult to trace. It is therefore reassuring that the supernova data indication of dark 
energy is confirmed by cosmic microwave background (CMB) anisotropy data. More pre- 
cisely, assuming the CDM model of structure formation and assuming that dark energy 
does not evol ye in time, CMB anisotropy data are consistent with negligible space curva- 
ture (see, e.g.JPodariu et ali hoOlbl; IPaee et alJliooS boran et al.1l2007l : bunklev et aDl2008 



Komatsu et al 



v see, e.g., 



200 81), so in con junction with low measured non-relativistic matter density 



Chen fc Ratrall2003bl ) CMB anisotropy data also demand dark energy. 



To convincingly remove degeneracies and understand and cancel the effect of unwanted 
systematic errors it's necessary to have many independent cosmological tests. Other re- 
cently discussed cosmological pr obes include the angu l ar size of radio sources and quasars as 



a function of redshift (see, e.g.. IChen fc Ratrall2003al ; iPodariu et al. 



Santos fc Lima 



2008 



2007 



Zh ang et al 



200 81), strong gravitationa 



2009 



lensing (see, e.g., Lee &: Ngl 



2003 



Zhu fc Serenol 2008 ), weak gravitational lensing ( 



Da 



2007 



v et al. 



2007 



Qguri et al. 



La Vacca fc Colombo! 120081 : iHoekstra fc Jainl I2OO8 



the Hubble parameter as a function of redshift (see, e.g., 



20081 : bev et alJliooi : iFernandez-Martinez fc Verde! 200 8]), galaxy cluster gas mass fra c tion 



versu s redshift data (see, e.g., 



Allen et al 



2004 



see, e.g. 



Pore et al. 



Schmidt 2008). measurements of 



Samushia fe Ratra 2006; Lin et al. 



Chen fc Ratra 



2004 



2008) , and the grow t h of large-scale struc ture (see, e.g.. lBasilakos et al. 
200a lMaininilEo08l : kbramo et al1l2008l ). 



Allen et al. 



2008 



2008: 



Sen 



Waizmann fc Bartelmann 



A measurement of the large-scale structure baryon acoustic oscillation (BAO) peak 
length scale provides another possible cosmological test. This was first measured a few years 
ago in the two-point correlation function of luminous red galaxies in the SPSS catalog at 
z = 0.35 (jEisenstein et al.ll2005l . also see iHutsil 120061 ) and from 2dFGRS da ta (ICole et al. 
2005). and then later in the joint SPSS and 2dF catalogs at z = 0.24 and 0.43 (jPercival et al.l 



2007a! ) H One attractive feature of the BAO measurements is that they do not explicitly 
depend on the value of Hubble constant. These BAO measurements wer e used in the o rigina l 
papers to constrain parameters of the ACPM and XCPM models. See llshida et al.l (120081 ). 



4 See Blake et al. ( 2007 ). Padmanabhan et al. ( 2007), an d Percival et al.l J2007bf) for more recent discussio n 
of the observational situation. See Blake k, Glazebrook ( 20031) . Linderj j2*003l ), Seo fc Eisenstein ( 2003 ). 



Hu fc Haimanl (|2003l ). and references therein, for discussions of the BAO peak length scale "standard ruler" 



cosmological test. 
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Lazkoz et al.l (120081 ) , ISantos fe Jesua (120081 ) , and ISamushia fe Ratral (j2008bl ) , and references 
therein, for recent discussions of BAO data constraints on these and other models. 



The iPercival et al.l (12007a ) measurement was made at only two redshifts and by itself 
does not provide a robust and reliable test of dark energy models, but in combination with 
other data it does provide useful constraints on cosmological parameters. However, a number 
of surveys are planned in next few years that will measure the BAO scale accurately and at 
a variety of redshifts up to z = 1.2. This upcoming BAO data, especially when combined 
with other data, will p rove very useful in tightly constraining dark energy parameters (see, 



e.g. 



Wang et all 120091 ). 



Recently IGaztanaga et al.l (j2008al ) argued that these measurements of the BAO scale 
were essentially measurements orthogonal to the line of sight and so statistically independent 
from a line of sight mea surement of the BAQ sca le, even if the same galaxy catalog is used 
for both measurements. IGaztanaga et al.l (j2008al ) used the SDSS data to compute the line 



of sight or radial BAO scale in redshift space for two ranges of redshift and IGaztanaga et al. 



( 12008bl ) showed the resulting constraints on the spatially-flat XCDM model. These constrains 
from the radial BAO scale data are quite similar to the constraints derived from earlier "non- 
radial" BAO scale measurements. In this paper we extend the analysis to the ACDM and 
0CDM models of dark energyjf] The data are consistent with spatially-flat ACDM. However, 
these current radial BAO measurements (like current non-radial BAO measurements) can not 
tightly constrain time-varying dark energy, although the situation is anticipated to improve 
in the next few years. 

We also derive constraints on th ese models from a combined analysis of the radial BAO 
data and the iKowalski et al.l (120081 ) Union type la supernova data. Since the radial BAO 
and type la supernova data constraints are almost orthogonal, the constraints from the 
combined data are significantly tighter than those from either individual data sets. The 
combined constraints favor a close to spatially-flat ACDM model (more so than the type la 
supernova data), but do not yet completely rule out time- varying dark energy. 

Our paper is organized as follows. In Sec. 2 we describe the radial BAO measurements 
we use. In Sec. 3 we explain how we derive constraints on different dark energy models from 
radial BAO and type la supernova data. We present and discuss our results in Sec. 4. 



IGaztanaga et all (|2008bl ) quote limits on the ACDM parameters but do not show constraint contours 
for ACDM. 
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2. Radial baryon acoustic scale 

In a spherically symmetric Universe the two-point correlation function is a function 
of two variables, £ = £(er, 7r), where a is the separation along the line of sight and rr is 
the separation on the sky. It can also be expressed as a function of absolute separation 
r = \/ a 2 + tt 2 and the cosine of the angle between the line of sight and the direction of 
separation, \l = n/r. The correlation function can then be decomposed into multipole 
moments 



6(r) = J i C(r,fi)Pi(fi)dn, (1) 

where Pi is the / th order Legendre polynomial. Multipole moments of different orders can 
be related to each other if one has a complete theory of linear and nonlinear evolution. 
Although high multipoles that describe the "shape" of baryon acoustic oscillation imprints 
on the matter distribution are very difficult to measure in practice, theoretically they are 
independent of the monopole and could provide additional structure formation tests. 

Initial work considered only the averaged over direction monopole part of the correlation 
function, 



£o(r) 



+i 



(2) 



and found a BAO peak signal at a comoving distance of r ~ 110/i _1 Mpc (h is the Hubble 
constant Hq in units of 100 kms _1 Mpc _1 ). This measurement was however mostly transverse 
to the lin e of sight direction 7r; the weight of separation along the line of sight contributes less 
then 1% (IGaztanaga et al.ll2008al ). Consequently, it is fair to assume that the radial baryon 
acoustic peak scale measurement in the line of sight direction from £(cr) is statistically 
independent from that measured from £(r) ~ £(7r), even if the same galaxy sample is used 
for both measurements. 



Gaztanaga et al.l (J2008aJ) recently used SDSS data to measure the radial baryon acoustic 
scale in two redshift ranges z ~ 0.15 — 0.30 with radial BAO peak scale Az = 0.0407± 0.0014 
and z ~ 0.40 — 0.47 with Az = 0.0442 ± 0.0016 (both one standard deviation errors). 
Theoretically the radial BAO peak scale is given by 



Az = H(z)r s (z d )/c 



(3) 
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where H(z) is the Hubble parameter at redshift z, r s (z<i) is the sound horizon size at the drag 
redshift z d , at which baryons were released from photons, and c is the speed of light. H(z) 
can be easily computed in a given cosmological model and depends on model parameters 
such as the non-relativistic matter density and the time dependence of dark energy 



To compute r s IGaztanaga et al.l (l2008bl ) use two different methods. One is to use the ra- 
tio l s between the distance to the last-scattering surface and r s measured by CMB anisotropy 
experiments and compute the sound horizon at photon decoupling from 

r.(s.) = ?r(1 + f^ ) . (4) 

Here z* is the redshift at the photon decoupling and (1a is the angular diameter distance. 
Alternatively, one can use priors on the fractional energy density parameters of baryonic mat- 
ter, f2b, nonrelativistic matter, fl m , and relativistic matter, Q T , from, e.g., CMB anisotropy 
measurements, and compute the sound horizon at the drag redshift from 

r s (z d ) = °-= r (Zd) da . (5) 

H ^TU Jo y/( a + 1.69f2 r /n m )(l + a0.75Q h /Q r ) 

Both options have similar drawbacks. One has to assume priors on "nuisance" param- 
eters like l s or various energy densities. CMB anisotropy measurements themselves have 
measurement errors that must be accounted for, otherwise the errors on the estimates of 
dark energy model parameters of interest will be underestimated. Also, the best fit values 
for nuisance parameters given by CMB anisotropy data are in general different for every cos- 
mological model and also depend on model parameter values. To be fully consistent when 
using priors one would have to reanalyze CMB experiments for each cosmological model (and 
model parameter value) instead of using a single set of values for l s , fl^, Q m , and Q r . 

At present, however, the BAO scale is measured only in two redshift ranges and does 
not provide very tight parameter constraints compared to other observational tests. Hence, 
as long as we are interested in preliminary constrai nts on dark ener g y from BAO scale 



measurements we may use the simplified approach of IGaztanaga et al.l (l2008bl ). keeping in 
mind that when more and better quality BAO scale measurements become available a more 
complete, careful, and time-consuming analysis will be warranted. 
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3. Computation 



We compare the two measured radial BAO peak scales (IGaztanaga et al.l l2008al ) with 
the predictions of three dark energy models. The models we consider are standard ACDM, 
the XCDM parametrization of the dark energy's equation of state, and the 0CDM model 
with an inverse power law potential energy density V(<ft) oc <p~ a . In these dark energy models, 
at late times, we can compute the redshift-dependent Hubble parameter from 



H(z) = tfoVMl + z f + + (1 - - Oa)(1 + z? (ACDM), (6) 

H{z) = H ^Q ni (l + zf + (1 - n m )(l + ;z)3(i+"*) (XCDM), (7) 

H{z) = #o^/On(l + z) 3 + fy(a,z) (0CDM). (8) 

In the 0CDM model the energy density of the scalar field has to be computed as a function 
of the redshift and a parameter by numerically solving the equations of motion^ 

In all three models the background evolution is described by two parameters. One is 
the nonrelativistic matter fractional energy density parameter Q m and the other one is a 
parameter p that characterizes the dark energy. For the ACDM model p is the cosmological 
constant fractional energy density parameter in the XCDM parametrization it is the 
equation of state parameter u x , and in the 0CDM model it is a which describes the steepness 
of the scalar field self-interaction potential. In this paper we consider only spatially-flat 
XCDM and 0CDM models, while in the ACDM case spatial curvature is allowed to vary, 
with the space curvature fractional energy density parameter i\ = 1 — Q m — Q^. The angular 
diameter distance is determined in terms of the Hubble parameter through 



d A (z) 
d A (z) 
d A (z) 



1 



■ sm 



^KH {1 + z) 

i r dz' 



dz' 
H(z' 



(l+z)J H(z') 
1 



V^kH {l + z] 



sinh 



V^Ho [ 
Jo 



dz' 
H(z' 



(fik > 0), (9) 
(O k = 0), (10) 
(fi k < 0). (11) 



We use both methods proposed in IGaztanaga et al.l (l2008bl ) to estimate the sound hori- 
zon at recombination. We compute theoretically predicted model values of the radial baryon 



3 For a discussion of the scalar field dynamics see, e.g., IPodariu fc Ratral |200d ) 



- 8- 



acoustic scale in redshift space, Az t h(^m,P, at redshift Z\ = 0.24 and z 2 = 0.43, where 
by v we denote the "nuisance" parameters. 

We assume that the measurements at z — 0.24 and 0.43 are independent and that the 
errors are Gaussianly distributed. To constrain model parameters we compute 

X 2 (fi m ,p, v) = ^2 (Az t h, t (fim,P, v) - Az ohStl ) 2 / a 2 AZti , (12) 

i=l,2 

where the two observed values, Az \, Sti ±(jAz,i, are listed above Eq. (3) and define a likelihood 
function 



£ B Ao(^m,P, v) cx exp(-x(fi m ,p, i/)/2). 



(13) 



We integrate over nuisance parameters to get a two dimensional likelihood for the cos- 
mological parameters of interest £bao(^iii,p) = / ^BAo(^m,P, v)L{y)dv. For the prior 
likelihood of nuisance parameters we use gaussian distribution functions with WMAP 5- year 
recommended means and variances l s = 302 ± 0.87, z s = 1090 ± 1 (with 40% positive cor- 
relation between t wo measurements), Q y,h 2 = 0.0227 ± 0.0066, Q m h 2 = 0.133 ± 0.0064, and 
Q r = 2.45 x 10 -5 ( iDunkley et al.ll2008l ). For each dark energy model we define the best- 
fit values of parameters as a pair (Q m , p) that maximizes the two-dimensional likelihood 
function £ B Ao(^m,£>)- The 1, 2, and 3a confidence level contours are defined as the sets 
of parameters (ft m , p) for which £ B Ao(^m,p) is less than its maximum value by multiplica- 
tive factors of exp(— 2.30/2), exp(— 6.18/2), and exp(— 11.83/2), respectively. In Figs. 1-3 
we show two-dimensional constraints on cosmological parameters from radial BAO data for 
each of three dark energy models and two sets of priors. 

Since radial BAO measurements alone, at the moment, can no t tight ly constrain cos- 
mological parameters we combined them with the iKowalski et al.l ( 120081 ) Union type la 
sup ernova data. I n our derivation of the type la supernova constraints we closely fol- 
low IKowalski et al.l (120081 ) . We define a likelihood function for cosmological parameters, 
£sN(^m>P> -^o)> f° r each dark energy model. We marginalize over H with a flat noninfor- 
mative prior in the range 40 < Ho/(km s _1 Mpc -1 ) < 100 to get a two-dimensional likelihood 
C^s(Q m ,p). The constraints on the 0CDM model from the supernova data are shown in Fig. 
4[] Since radial BAO and supernova measurements are independent, we multiply the indi- 
vidual likelihoods to get a joint likelihood function C(Q m ,p) = C B x (Q m ,p)C S ' N (Q m ,p). We 



7 For the constrain ts on ACDM and XCDM models from the supernova "Union" data set see Fig. 12 in 
Kowalski etall |2008l ). 
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show constraints derived from the joint likelihood function in Figs. 5-7. 

From the joint two-dimensional likelihood function we define the marginal one-dimensional 
likelihood functions for parameters Q m and p (assuming uniform priors on both parameters 
in the ranges Q m G (0.0, 1.0), fi A G (0.0, 1.0), u x G (-2.0, 0.0), and a G (0.0, 5.0)) through 



C(tt m ) = j C(Q m ,p)dp, (14) 

c{p) = J c(n m , P )dn m . (15) 

In Fig. 8 we show one-dimensional likelihood functions for fl m and p for the three dark 
energy models. From these one-dimensional likelihood functions, in all three models, we 
define the best-fit value of parameters as the values that maximize the likelihood, with 
highest posterior density la confidence level intervals as the value of parameters for which 
J ' la . dx C(x) j fy x dx C(x) = 0.68 and with C(x) higher everywhere inside the interval than 
outside. The values of one- dimensional best-fit parameters and corresponding la confidence 
level intervals are listed in Table 1. 



4. Results and discussion 



Figure 1 shows the lGaztanaga et al.l (J2008a]) radial BAO scale constraints on the ACDM 
model. When we use the WMAP measured value of the ratio l s (thick lines) the model is 
constrained to be very close to the spatially-flat case. When we use the WMAP measured 
value of the sound horizon at recombination r s (thin lines) the constraints are weaker, the 
spatial curvature is not well constrained, and the nonrelativistic matter Q m has to be less than 
0.45 at about 3a. The c ontour s in Fig. 1 are in reasonable accord with those shown in Fig. 
2 of ISamushia &: Ratral (I2008bf) wh ich were derived using the " non-radial" BAO peak scales 
measured by lPercival et al.l (j2007al ) and lEisenstein et al.l (120051 ) . When the radial BAO data 
is combined with supernova data the resulting constraints are significantly stronger, see Fig. 
5. The nonrelativistic matter density parameter is now constrained to be 0.15 < Q m < 0.35 
at about 3a confidence level while the cosmological constant density parameter lies in the 
0.45 < Q\ < 0.9 range. The best-fit values are close to the spatially-flat model. The 
constraints computed from the joint one- dimensional likelihoods shown in Figs. 8a and 8b 
are listed in Table 1. They are much more restrictive than those derived from the individual 
data sets, because those constrain combinations of dark energy parameters that are in some 
sense "orthogonal" in parameter space. 
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For the XCDM parametrization (for which we consider only spatially-flat models) the 
confidence level contours derived from radial BAO measurements are broad and a range of 
Q m and uj x values ar e acceptable. The constr aints are shown in Fig. 2. They are similar to 
the results shown in iGaztanaga et al.l (l2008bl . Fig. 1). When we use the WMAP measured 
value of the ratio l s (thick lines) the nonrelativistic matter Q m has to be less than 0.4 at 
about 3a. These confidence level contours are in reasonable accord with the ones shown 



m iSamushia fc Ratral (I2008bl. Fig. 3) which were d erived using the non-radial BAO data 
of Percival et ali i 2007a ) and Eisenstein et al. ( 2005 ). However, compared to the non-radial 
BAO scale measurements the radial BAO scale measurements better constrain u x from below 
and tend to favor higher values of it. The joint constraints from radial BAO and supernova 
data are shown in Fig. 6. The joint likelihood constrains the equation of state parameter to 
be —0.7 < lj x < —1.2 at about 3a confidence, while 0.2 < Q m < 0.35. For XCDM the joint 
constraints do not depend much on the method used for the analysis of the radial BAO data, 
unlike the ACDM and 0CDM cases. One-dimensional joint likelihood functions for fl m and 
u x are shown in Figs. 8c and 8d and the corresponding constraints are given in Table 1. 

The spatially-flat 0CDM model confidence level contours are shown in Fig. 3. Here, the 
radial BAO measurements constrain f2 m to be between 0.15 and 0.4 at about 3a, but the 
a parameter is not constrained well and large values of a (relatively rapidly evolving dark 
energy) are not rule d out, although the l ikeliho od peaks at a = 0. These results are similar to 
the ones derived in lSamushia fc Ratral (I2008bl. Fig. 4) using t he no n-radial BAO peak scale 



measurements of Percival et al. 



(j2007al ) and lEisenstein et al. J2005h . Figure 6 shows type la 
supernova constraints on 0CDM. From the Union data set alone, the nonrelativistic matter 
energy density parameter Q m < 0.4, and if fl m > 0.1, a < 4.3 at about 3a confidence. 
The se contours are si milar to and a little more const raining than the on es derived using 
the iRiess et al.l (120041 ) type la supernova Gold data (I Wilson et al.l 120061 . Fig. 1). When 
we combine the radial BAO measurements with the Union supernova data, Fig. 7, the a 
parameter is constrained to be less than 1.5 at 3a while 0.2 < Q m < 0.35. One- dimensional 
joint likelihood functions for fl m and a are shown in Figs. 8e and 8f and the corresponding 
constraints are listed in Table 1. For both radial BAO analysis methods fl m is relatively 
well constrained and the likelihood for a peaks at low values of a, close to the corresponding 
spatially-flat ACDM values. 



From the analysis presented here and in ISamushia fc Ratral (j2008bl ). it is clear that 
neither the current radial or non-radial BAO peak scale measurements by themselves are 
able to tightly constrain cosmological parameters. However, given the shape of the resulting 
confidence contours, typically an elongated ellipse-like shape, current BAO data can provide 
very useful constraints when combined with other data like type la supernova measure- 
ments. As another example, the constraints derived from radial and non-radial BAO peak 
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meas urements complement the constraints derive d from galaxy cluster gas mass fraction 
data (lAUen et al.l 120081 ; ISamushia fc Ratral l2008af) and when used join tly they can tightly 
constrain parameters of dark energy models (ISamushia fc Ratral l2008bl ). It is interesting to 
note that for the ACDM and XCDM models the best fit values from the BAO and g alaxy 



cluster gas mass fraction data sets are about 3a apart (see ISamushia fc Ratral l2008bl . Figs. 
2-3), but for the two data sets together the minimum x 2 is acceptable. This is, most proba- 
bly, due to unknown systematic errors in one (or both) of the data sets and deserves further 
investigation. 

It is also of interest to compare the constraints shown in Figs. 1-3 to those derived using 
other data. The radial BAO constraints shown in Figs. 1 -3 are more restrictiv e than the 
ones derived from Hubble parameter versus redshift data (ISamushia et al.l 120071 ) and from 



gravitational lensing data (jChae et al.l 12004 



derived from earlier type la supernova data (IRiess et al.l 120041 ) by lWilson et al.l (120061 ). 



However, they ar e comparable to the ones 



Current BAO peak scale data is sparse and can not tightly constrain dark energy param- 
eters or differentiate betw een different dark e nergy models. In nex t few years a number o f 
planned surveys like PAU faenitez et al.ll2009h . BOSS, and WiggleZ Jdazebrook et alil2007f ) 
should measure the BAO peak scale with 1% accuracy possibly up to a redshift of z = 1.2. 
These future BAO measurements should prove to be of great significance in discriminating 
between different currently viable dark energy models and constraining cosmological param- 
eters. 



We acknowledge useful discussions with E. Gaztanaga and support from DOE grant 
DE-FG03-99EP41093 and INTAS grant 061000017-9258. 
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Fig. 1. — 1,2, and 3cr confidence level contours for the ACDM model from radial BAO mea- 
surements. Thick lines are derived using the WMAP measured ratio l s and the cross shows 
best-fit model parameters, Q m = 0.21 and = 0.80, derived from the two-dimensional like- 
lihood function. The thin lines are derived using the WMAP value for the sound horizon at 
recombination r s and the circle shows best-fit model parameters, Q m = 0.23 and = 0.80, 
derived from the two-dimensional likelihood function. The dotted line demarcates spatially- 
flat ACDM models. 
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Fig. 2. — 1, 2, and 3a confidence level contours for the XCDM parametrization from radial 
BAO measurements. Thick lines are derived using the WMAP measured ratio l s and the 
cross shows best-fit model parameters, Q m = 0.22 and u x = —1.03, derived from the two- 
dimensional likelihood function. The thin lines are derived using the WMAP value for the 
sound horizon at recombination r s and the circle shows best-fit model parameters, Q m = 
0.22 and f2 x = —1.03, derived from the two-dimensional likelihood function. The dotted 
horizontal line demarcates spatially-flat ACDM models. 
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Fig. 3. — 1, 2, and 3a confidence level contours for the 0CDM model from radial BAO 
measurements. Thick lines are derived using the WMAP measured ratio l s and the cross 
shows best-fit model parameters, Q m = 0.23 and a = 0.0, derived from the two-dimensional 
likelihood function. The thin lines are derived using the WMAP value for the sound horizon 
at recombination r s and the circle shows best-fit model parameters, Q m = 0.24 and a = 
0.0, derived from the two-dimensional likelihood function. The horizontal axis with a = 
corresponds to spatially-flat ACDM models. 



-19- 




Fig. 4. — 1, 2, and 3a confidence level contours for the 0CDM model from the supernova 
data. The cross shows best-fit model parameters, Q m = 0.29 and a = 0.0, derived from 
the two-dimensional likelihood function. The horizontal axis with a = corresponds to 
spatially-flat ACDM models. 
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Fig. 5. — 1, 2, and 3er confidence level contours for the ACDM model derived using the 
joint likelihood from radial BAO and type la supernova measurements. Thick lines are 
derived using the WMAP measured ratio l s and the cross shows best-fit model parameters, 
Q m = 0.27 and Q\ = 0.70, derived from the two-dimensional likelihood function. The thin 
lines are derived using the WMAP value for the sound horizon at recombination r s and the 
circle shows best-fit model parameters, Q m = 0.25 and Q\ = 0.68, derived from the two- 
dimensional likelihood function. The dotted line demarcates spatially-flat ACDM models. 
The fi m and ranges shown here are smaller than those shown in Fig. I. 
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Fig. 6. — 1, 2, and 3a confidence level contours for the XCDM parametrization derived using 
the joint likelihood from radial BAO and type la supernova measurements. Thick lines are 
derived using the WMAP measured ratio l s and the cross shows best-fit model parameters, 
Q m = 0.24 and u x = —0.91, derived from the two-dimensional likelihood function. The thin 
lines are derived using the WMAP value for the sound horizon at recombination r s and the 
circle shows best-fit model parameters, Q m = 0.25 and fl x = —0.93, derived from the two- 
dimensional likelihood function. The dotted horizontal line demarcates spatially-flat ACDM 
models. The f2 m and u x ranges shown here are smaller than those shown in Fig. 2. 
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Fig. 7. — 1, 2, and 3<x confidence level contours for the 0CDM model derived using the 
joint likelihood from radial BAO and type la supernova measurements. Thick lines are 
derived using the WMAP measured ratio l s and the cross shows best-fit model parameters, 
Q m = 0.24 and a = 0.27, derived from the two-dimensional likelihood function. The thin 
lines are derived using the WMAP value for the sound horizon at recombination r s and the 
circle shows best-fit model parameters, Q m = 0.25 and a = 0.20, derived from the two- 
dimensional likelihood function. The horizontal axis with a = corresponds to spatially-flat 
ACDM models. The tt m and a ranges shown here are smaller than those shown in Figs. 3 
and 4. 
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a) ACDM b) ACDM 




Fig. 8. — One-dimensional likelihood functions for cosmo logical parameters for the three 
models derived using the joint likelihood from radial BAO and type la supernova data. 
Thick lines show results derived using the WMAP measured ratio l s while thin lines are 
derived using the WMAP value for the sound horizon at recombination r s . In each of the 
six subplots total likelihoods are normalized to one. 
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Table 1. Best fit values and la range of Q m and p. a 



model 




Q m la range 




P 


p la range 


ACDM (prior l s ) 


0.27 


0.24-0.31 




= 0.70 


0.64 < tt A < 0.75 


ACDM (prior r s ) 


0.25 


0.24-0.32 


n A 


= 0.68 


0.62 < tt A < 0.77 


XCDM (prior l s ) 


0.24 


0.23-0.27 


u x = 


= -0.91 


-1.00 < uj x < -0.88 


XCDM (prior r s ) 


0.25 


0.23-0.28 


W s = 


= -0.93 


-1.03 < uj x < -0.85 


0CDM (prior l s ) 


0.24 


0.22 - 0.26 


a - 


= 0.27 


0.02 < a < 0.54 


0CDM (prior r s ) 


0.25 


0.24-0.27 


a - 


= 0.20 


a < 0.43 



a Best-fit values of cosmological parameters Q m and p, and la confidence level intervals, 
from one-dimensional likelihood functions derived from a joint analysis of radial BAO and 
type la supernova data. Entries labeled as "prior l s " are computed using the WMAP mea- 
sured ratio l s (and correspond to the thick lines in the figures). Entries labeled as "prior r s " 
are computed using the WMAP measured value for the sound horizon at recombination r s 
(and correspond to the thin lines in the figures). 



